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ABSTRACT 



Recently, Buta et al. (2009, AJ, 137, 4487) examined the question "Do Bars 
Drive Spiral Density Waves?" , an idea sup ported by theoretic al studies and also 



200J). They estimated 



from a preliminary observational analysis ( iBlock et al. 
maximum bar strengths Qb, maximum spiral strengths Qs-, and maximum m=2 
arm contrasts for 23 galaxies with deep AAT fr^-band images. These were 
combined with previously published Qh and Qs val ues for 147 ga l axies from the 



Block et al. 



f l2004h . Weak 



OSUBSGS sample and with the 12 galaxies from 
correlation between Qj, and Qs was confirmed for the combined sample, whereas 
the AAT subset alone showed no significant correlations between Qh and Qs, nor 



between and A2s- A similar negative result was obtained in 



Durbala et al. 



( I2OO9I ) for 46 galaxies. Based on these studies, the answer to the above question 
remains uncertain. 

Here we use a novel approach, and show that although the correlation be- 
tween the maximum bar and spiral parameters is weak, these parameters do 
correlate when compared locally. For the OSUBSGS sample a statistically sig- 
nificant correlation is found between the local spiral amplitude, and the forcing 
due to the bar's potential at the same distance, out to ^1.6 bar radii (the typical 
bar perturbation is then of the order of a few percent). Also for the sample of 



23 AAT galaxies of 



Buta et al. 



(I2OO9I ) we find a significant correlation between 
local parameters out to bar radii. Our new results confirm that, at least in 
a statistical sense, bars do indeed drive spiral density waves. 



Subject headings: galaxies: spiral — galaxies: kinematics and dynamics — galaxies: 
structure 
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Introduction 



The response of an outer disk to a rotating bar is intima tely relatec 



of long lasting spiral arms in galaxies. Optical photometry (iSchweizer 



to the maintenance 



19761) established 



that besides the gas and young stars, the spi rals are present also in the old population 



This was confirmed by n ear-infrared surveys (lEskridge et al. 



for grand-design spirals (IKnapen fc Beckman 



be induced via disk instabilities fiBertin et al. 



galaxy interactio ns (iToomre fc Toomre 



galaxy rotations ( ISellwood fc Carlberg 



cycle, e.g. due to the swing amplification f lToomre 



2OO2I ). and is particularly true 



9961) . Short-hved stellar density w aves can 



1977 



Goldreich fc Tremaine 



1978h . or by 



1972), but such transient patterns fade after ~10 



1984). unless they are maintained by some feedback 



19811 ). On the other hand, spiral arms 



are e xcited by a growing bar, as dem onstrated by the ver y first N-body simulations (IHohl 



I971I ) and by analytical calculations (lAthanassoula 



19801 ). 



Support for the bar/spiral connection is provided by the examples where prominent 
spirals extend from the ends of the bar (see e.g. NGC 1300, p. 525 in Binney & Tremaine 
2008, or NGC 986 in Buta et al. 2 010). Also, grand-design spirals are more frequent in 



barred than in non-barred galaxies ( lElmegreen fc Elmegreen 



19821 ). Nevertheless, a direct 



connection between bars and spirals has been difficult to prove observationally. 



The possible driving of spirals by bars was addressed by iBlock et al. 



(120041 ) ■ who 



compared the maximum of bar-related torque strength Qb (the maximum of tangential 
force amplitude normalized by mean radial force) with the maximum associated with the 
spirals ,Qs, after separating the bar and spiral components based on their Fourier density 
amplitude profiles. Near-IR observations were used, to show the effect of the bar on the 
surrounding mass, rather than on the gas or the form ation of young stars. Based on 12 



galaxies observed in the i^s-band. 



Block et al. 



( I2OO4J ) found a strikingly clear correlation 



between Qb and Qs- However, the correlation is less clear in later studies using larger 
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samples and deeper images. In iButa et al. I (120051 ) we analyzed the if-band images from 
the Ohio State University Bright Spiral Galaxy Survey (OSUBSGS, Eskridge et al. 2002) 
for 147 galaxies, and were able to "weakly verify a possible correlation between Qg and 
Qb" ■ More recently, in 



Buta et al. 



( l2009l ) we analyzed deep Kg Anglo-Australian Telescope 
(AAT) observations for 23 galaxies, with no statistically significant correlation between 
Qh and Qs, nor betw een Qh and (the m aximum of m=2 Fourier density amplitude 



of spirals). Similarly, 
z-band data for 46 isolated 
strengths wa s seen also by 



Durbala et al. 



( 120091 ) found no correlations when analyzing Sloan 



jarred galaxies; a la ck of correlation between bar and spiral arm 



Buta et al. 



Seiger et al. 



(|2003[ ) who analysed 41 galaxies. Nevertheless, 



( I2OO9I ) showed th at the correla t ion is pres ent when combinin g the AAT data 



with the previous data sets of 



Block et al. 



(120041 ) and 



Buta et al. 



(l2005h . 



In this p aper the bar/ s piral conne ction is re- investigated using the same samples which 



were used by 



Buta et al. 



( 12005 



20091 ): the OSUBSGS sample containing ~100 barred 



galaxies, and our AAT sample of 23 barred galaxies. A novel approach is used: instead 
of comparing the maximum bar strength with the maximum spiral density amplitude, we 
compare the locally measured bar forcing and spiral amplitude as a function of distance. 
The locally felt forcing due to bar is a more important parameter than the maximum 
forcing, since the Qt, is typically attained well inside the spiral structure. Also, we examine 
the spiral density rather than force amplitude, since the former measures more directly the 
possible response to bar forcing. Using our approach, a statistically significant correlation is 
demonstrated to exist between the bar forcing and the spiral amplitude, up to a considerable 
distance beyond the end of the bar. 
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2. Calculation of bar forcing and spiral amplitudes 

We calculate the amplitude of the bar tangential forcing as a function of distance 
from the galaxy center, Qtari^), and compare it to the m=2 surface brightness amplitude 
of the spirals, A2{r), at the same radial distances, normalized to bar radius, r/Ri,ar- The 
Spearman rank correlation coefficient between Qbar and A2 and its significance level are 
then measured as a function of r/Rbar- A non-parametric test is used in order to avoid 
making assumptions about the distributions of the compared quantities. The distance up 
to which a significant correlation is found is regarded as a statistical estimate of the region 
inside which bars are able to drive spiral structure. 

We use near-IR images to evaluate the galaxy potential, and derive the tangential force 
amplitude at each distance, normalized to the azimuthally averaged radial force 



Qbar{r) 



max(|FT(r, (/?)!) 



< |FR(r,(^)| > 

Several assumptions are made: (1) the mass-to-luminosity ratio (M/L) is constant, (2) the 
vertical profile of the disk and the bar is approximated with an exponentia l function, and (3 ) 



the v ertical scale height, hz, scales with the galaxy size as hz = O- lRr^on (ISpeltincx et a 



20081 ). where Rk2o is the rriK = 20 isophotal radiu s from 2MASS (ISkrutskie et al. 



The calculations are made with polar integration (iSalo et al. 



199f 



2006) 



Laurikainen fc Salo 



2OO2I ). based on azimuthal Fourier decomposition of the deprojected image 



J(r,(^) = Ao{r) 



1 + ^ Am{r) cos (m - ipm{r)]) 



m=l 



which also provides the m = 2 Fourier density profile used to characterize the spiral 
amplitudes (note the normalization). The gravity is calculated separately for each m 
component (using m = 0,2,4,6,8, 10) and added together. Compared to direct Cartesian 
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integration, the polar method suppresses effect ively the spurious force maxima which 



otherwise could arise in the noisy outer disks ( ISalo et al. 



20041 ) ■ To account for the 



different 3D-distribution of the bulge light, mult i- component decompositions are used 



( iLaurikainen et al. 



20051 ). reducing the artificial tangential force amplitudes arising from the 
bulge deprojection stretch. The polar method makes it also easy to exclude the contribution 
of spirals on the calculation of Qban by setting the m > Fourier density amplitudes to 
zero beyond a certain cutting distance, Rcut, representing the end of the bar. 

Our method is illustrated in Figure 1, where the radial profiles of the m=2 and m=4 
Fourier density amplitudes, and the Qbar ratio are shown for NGC 1566. In this particular 
example the A2 and have two well-separated local maxima associated with the bar and 
spiral, but the distinction is not always clear. The Qhar profile is constructed both for the 
total (bar+spiral) force field, and without the contribution of the spiral arms. The bar-only 



profile is obtained using a cutting distance Rcut = R hnr in the force calcu . 



But a et al. 



( 12003 



ation. Our method 



2005 



20091 ). who 



of isolating the bar forcing is different from that of 
extrapolated the bar density into the spiral region based on Gaussian fits to Fourier density 
profiles. Here we assume that the bar and the spiral dominate their radial domains with no 
significant overlap. 



3. Bar driven spiral structure in the OSUBSGS and AAT samples 



The OSUBSGS is a magnitude-limited sample (m^ < 12.0 mag) of galaxies with 
Hubble types < T < 9. The //-band images typicall y reach 20 mag/arc s ec^ in depth. 



Laurikainen et al. 



mm . and are 



Our bar identifications and lengths are from Table 3 in 
based on the Fourier amplitude and phase profiles ('Fourier bars'; NGC 2207 was omitted, 
leaving 103 galaxies for analysis). 
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A statistically significant correlation is found between the amplitude A2, and the bar 
forcing, Qbar, when these parameters are compared at the same radial distances. See Fig. 2 
where the different panels represent examples of measurements at successively larger radial 
distances with respect to the bar. The bar force is cut using Rcut = Rbar, which means 
that we are conservative in eliminating any contamination by the spiral arms themselves 
in the forcing. The correlation is very strong just beyond the bar, and stays statistically 
significant until ^l.GRbar (rank correlation coefficient = 0.25, significance p = 0.008). 
The correlation is similar for long {Rbar/hR > 1) and short {Rbar/hR < 1) bars, when 
normalized to the disk scale length. The range of the significant correlation is similar also 
for early (T < 3) and late type (T > 4) spirals (not shown in the plots). 

Figure 3 collects the correlation coefficients between A2 and Qbar at different distances 
^/Rbar- Note that in the bar region, where A2 represents the bar itself, the correlation is 
strong as expected. Outside the bar A2 arises from spiral structures (the forcing is still 
dominated by the bar). The radial trend depends only slightly on the adopted cutting 
distance for the bar (using reasonable values Rcut/ Rbar = 1-0 — 1.2). The correlation is also 
insensitive to the exact assumptions made about vertical structure: an uncertainty by a 
factor of two in hz corresponds to ~ 20% uncertainty in bar strength. Monte Carlo trials 
show that even 20% random errors in Qbar (and A2) affect only marginally the significance 
of the correlation (for 10, 000 trials the median p = 0.023 for r = l.QRbar)- 

We made a simil ar analysis for the AAT sample of 23 barred galaxies, using the data 



and bar lengths from 



But a et al. 



( 120091 ). The i^'g-band images typically reach a surface 
brightness level of 22-24 mag/arcsec^. The sample covers a wide range of bar strengths, 
extending to strongly barred galaxies, for which ca se a stronge s t corr elation between 



Qbar and A2 is expected. However, as discussed in 



Buta et al. 



(l2009[ l. no statistically 



significant correlation is obtained between the maximum of Qbar and the maximum of 
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A2 (Fig. 4a). The result is similar for galaxies having the maximum spiral amplitude 
nearer (_R2s < l-S-R^ar) or further (i?2s > ^-^Rbar) outside the bar. However, a statistically 
significant correlation is found if the local bar forcing at the location of the maximum 
spiral amplitude is examined (Fig. 4b). The correlation is particularly clear if we limit to 
cases where i?2s < l-S-R^ar- This is in accordance with our result for the OSUBSGS where 
a statistical dependence is present but discernible only up to certain distance beyond the 
bar end. Indeed, if we repeat the analysis we made for the OSUBSGS sample (Fig. 4c), a 
statistically significant correlation is found between Qbar{r) and A2{r) up to 1.4 Rtar- The 
somewhat smaller range in the AAT analysis is probably due to the smaller sample size in 
comparison to the OSUBSGS sample. 



4. Discussion 



4.1. Comparison to previous studies of bar/spiral correlation 



Our analysis for OSUBSGS has indicated a significant correlation between the local 
tangential bar force Qbar{r) and the local surface brightness Fourier amplitude A2{r), up to 
l.GRbar- How does this compare with the previo us analyses of Qb y s. Qs, some of which 



But a et al. 



(l2005[ l we did not explicitly 



indicated a correlation, while others did not? In 
state the correlation coefficient, but using the data tabulated in the paper one finds a 
Pearson linear correlation coefficient r = 0.35. For a sample of = 146 galaxies this implies 
a very significant correlation {p = 0.8 ■ 10~^). Though different quantities are compared, 
this statistically significant correlation between maximum values agrees with our present 
analysis of the correlation between local quantities. 



Li 

and in 



Durbala et al. 



cewise, the negative r esults in iButa et al. I (l2009[ l for the AAT sample alone (N=23) 



( I2OO9I ) for the Sloan sample (N=46), are accounted for by the smaller 
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number of galaxies. The standard procedure for testing the significance of a positive hnear 
correlation (see e.g. Wall & Jenkins 2003, Sect. 4.2.2) implies that in order to accept the 
correlation (at a significance level p < 0.01), the sample correlation coefficients must be 
r > 0.48 and r > 0.34, for = 23 and = 46, respectively. Assume for a while that a 
linear correlation between Qb and Qs exists, and that these quantities are drawn from a 
bivariate-Gaussian distribution with an actual correlation coefficient p = 0.35 (equal to the 
sample correlation coefficient in Buta et al. 2005). We may then ask what the odds are 
of detecting this correlation, i.e. of observing a sufficiently high sample correlation when 
drawing a random sample with dif ferent N's. Applying t he Fisher probability distribution 



for the sample r with a known p (jWall fc Jenkins 



m i 



2003!) implies that for A^ = 23 there 



is only a 25% chance of detecting the correlation, even for A^ = 46 the change is only 
about 50%. If we use the rank correlation coefficient instead of the linear correlation 
coefficient, the chances are reduced to 15% and 38%, respectively (obtained by Monte Carlo 
trial estimates). Thus the negative results for these small samples do not rule out true 
correlations. 

The current method seems to be capable of exposing the correlation even for fairly 
small samples, indicating the advantage of comparing the local quantities. In particular, 
close to the bar the spiral amplitudes are strongly correlated with the bar forcing. This 
ra dial dependence of t he correlation probably explains the very strong correlation found 



m 



Block et al. 



( 120041 ). The fact that their small sample (N=12) showed the strongest 
correlation (r=0.86) is likely to result from the way the galaxies were selected, containing 
many examples where the spirals are strong right at the ends of the bar. Indeed, according 
to the tabulated values in their paper, the mean ratio between the radial locations of the 
spiral and bar ma ximum forces was < Rg/Rb >= 2.5, compared to < Rg/Rb >= 4.6 in 



Buta et al. 



( 120051 ) ■ which was based on a magnitude limited sample of spirals, with many 



different types of bar/spirals. 
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4.2. Physical mechanisms 



Our observational analysis indicates a clear statistical relation between bars and spirals. 
At least the following mechanisms might account for this: 

1) Spiral structure is a direct response to bar forcing and /or the spirals represent a 



continuation of the density wave associated with the bar ( iLindblad 



1960 



Toomre 



1969 



Bertin fc Lin 



19961 ) ■ Recently, spiral arms have also been int erpreted as manifold orbits 



emanating from unstable L agrangian points near the bar ends (iRomero-Gomez et al 



Athanassoula et al. 



2006 



20091 ) 



2) Spirals are 



Masset fc Tagger 



coup led to the bar via non-linear resonance coupling flTagger et al. 



19971 ). Such couplings are seen in N-body simulations (IRautiainen fc Salo 



1987 



19991 ). but it is uncertain how frequent such cases really are. 



The first explanation is likely to apply to the strong correlation just outside the bar. 
In this case the spirals are expected to share the constant (or slowly-evolving) pattern 
speed of the bar. Similarly, in the case of non-linear coupling, though the spiral pattern 
speed is slower than that of the bar, it still represents a steady long-lived pattern. On 
the other hand, in the region where no correlation is present, the spirals are independent 
stru ctures, representing either a long-lasting mode with a slower pattern speed than the 



bar (ISellwood fc Sparke 



of propagation speeds ( ISellwood fc Kahn 



1988r). or are just short-lived transient wave pack ets with a range 



1991 



Salo &: Laurikainen 



200' 



Note that some correlation may be present even if bars and spirals are independent, 
since both types of structures are favored by gravitationally more reactive disks. 
Nevertheless in this case a correlation between maximum values would also be expected. 



Sellwood 



(I2OO8I ) discusses the importance of distinguishing between long-lived spiral 



modes and transient waves. This stems from the fact that the latter are much more efficient 
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in driving angular momentum transport in the disk. In case of steady patterns the angular 
momentum exchange with stars is limited to resonances, whereas in the case of transient 
patterns this may occur over a large range of radii. Multiple transient patterns also lead to 
radial mixing of stars, as well as secular heating of the disk. 

In the current samples no significant correlation is seen beyond ~ 1.6Rf,ar, at which 
distance the typical force amplitudes associated with the bar fall below a few percent level. 
This provides an observational lower estimate for the extent of bar-driving. Namely, in the 
outer disk the determination of the spiral amplitude is more prone to uncertainties due 
to image noise and background subtraction, diluting any correlation that may be present. 
Also, we are fairly conservative in cutting the bar contribution close to Rbar, since in the 
case where the spiral arms are a part of the mode associated with the bar, the local spiral 
forcing would also contribute to maintaining the pattern. 

Although the correlation we find is strictly statistical, it is interesting to consider 
how it applies to the individual example of NGC 1566, shown in Fig. 1. For t his galaxy 



Rhnr = 40", and the A2 maximum related to the strong spirals is at = 68" (iButa et al. 



20091 ). though there are additional spiral arms even beyond 100". The strong spirals, with 
a maximum at 1.7 bar lengths, fall marginally in the region where bar-driving is expected. 
Continuing the same interpretation, the spiral beyond ~ 100" would be an independent 
pattern. Indeed, NGC 156 6 has been qu oted as an example of a galaxy with at least two 



separate spiral structures ( iBosma 



19921). 



5. Conclusions 



A connection between bar forcing and spiral density amplitudes was investigated 
for two near-infrared galaxy samples: 103 barred galaxies from the magnitude-limited 
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OSUBSGS survey, and 23 barred galaxies in our AAT Survey. The main results are: 

(1) For both samples a statistically significant correlation is found between the local 
tangential bar forcing, Qbar{f), o-nd the local spiral amplitude, A2{r), up to a radial distance 

Ofr r^l.SRbar- 

The correlation suggests that, at least in a statistical sense, the stellar spirals of the disk 
are not transient features but rather represent a continuation of the bar mode itself, or 
are driven by the bar through some mechanism. Further out, the spirals may be either 
independent modes or transient wave packets. 

(2) The obtained range of the correlation is similar for early and late-type spirals 

(0 <T <3 and 4: < T < 9) , and also for small and large bars (Rbar/W smaller /larger than 
unity). 

This does not favor the idea that only certain types of bars could drive spiral structure, or 
that the forcing on the stellar component requires the presence of significant gas component. 
Nevertheless, the current samples are small, and a larger number of galaxies is needed to 
draw reliable conclusions about morphological type dependencies, or to probe whether the 
statistical correlation extends to even larger distances beyond the bar. In this respect the 
forthcoming S^G survey (The Spitzer Survey of Stellar Structure in Galaxies; Sheth 2009, 
Sheth et al. 2010) will be extremely useful, providing unprecedentedly deep 3.6 and 4.5 
micron observations for nearly 2300 nearby galaxies. 
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Fig. 1. — Barred galaxy NGC 1566, with strong outer spirals. In a) the normalized A2 and 
A4 Fourier amplitudes are shown as a function of radius. The two local maxima are associated 
with the bar and the spiral: the vertical lines indicate the bar length Rbar, and the distance 
of maximum A2 of the spirals {R2s)- b) The calculated tangential forcing Qtar = FT/FR, 
with and without the contribution from spiral arms: in the latter case the m > Fourier 
components hav e been set to zero for r > Rcut, here using Rcut = Rbar- c) The deprojected 



Kg-hand image (IButa et al. 



20091 ). together with circles of radii Rbar and i?2s- 
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Fig. 2. — Relation between local bar forcing and local spiral amplitude. Four distances, 
normalized to the bar length {r / R^ar = 1.1 — 1.7) are compared, for the 103 barred OSUBSGS 
galaxies. In the calculation of bar forcing, the m > density Fourier amplitudes are set 
to zero beyond Rout — Rbar- The p values indicate the significance of the Spearman rank 
correlation coefficient r sample (the probability of having > r sample is P, under the hypothesis 
that the variables are independent). In the case where p < 0.01, the best-fit linear relation is 
also indicated. Open and filled circles denote short [Rbar/hu < 1) and long bars {Rbar/hR > 
1), showing no difference. 
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OSUBSGS: #barred gal =1 03 




Fig. 3. — Rank correlation coefficient between the local bar forcing Qbar{r) and the local 
amplitude A2{r), as a function of distance r/Rhar- Note that at r/Rhar < 1 the parameter 
A2 describes the bar density contrast whereas outside the bar it describes the spiral arms. 
Symbols indicate statistically significant correlation, obtained for r < l.BRbar- Note that the 
result is independent of the exact manner in which the spiral contribution has been eliminated 
from the bar forcing: compare the solid {Rcut = Rbar) and dashed {Rcut = ^■'^Rbar) curves. 
Also shown is the case with no cutting (dotted line): this is just for comparison, since the 
calculated forcing is then strongly affected by the spirals themselves. 
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Fig. 4. — Re-analysis of the AAT 



spiral components in iButa et al. 



samp 



e of 23 barred galaxies. The separation of bar and 
(120091 ) yielded the maximum of bar forcing Q^ar and the 
maximum of spiral amplitude A2, the latter attained at distance i?2s- Frame a) indicates no 
correlation between these quantities (same data as in Fig. 26a in Buta et al.). The filled and 
open symbols distinguish the galaxies with R2s/Rbar < 1-5 and R2s/ Rbar > 1-5, respectively, 
b) However, a correlation is present (with a significance p = 0.007) if we compare the bar 
forcing at the location of the maximum spiral amplitude; this is particularly clear for galaxies 
with R2s/ Rbar < 1-5 (filled symbols), c) The connection of bars and spirals is even clearer if 
the local bar forcing and local spiral amplitude are compared as a function of distance: the 
figure shows the rank correlation coefficient vs. distance, constructed in a same manner as 
in Figs. 2 and 3. 



